Effect of Disorder Strength on Optimal Paths 
in Complex Networks 
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We study the transition between the strong and weak disorder regimes in the scaling properties 
of the average optimal path £ opt in a disordered Erdos-Renyi (ER) random network and scale-free 
(SF) network. Each link i is associated with a weight n = exp(ar;), where n is a random number 
taken from a uniform distribution between and 1 and the parameter a controls the strength of 
the disorder. We find that for any finite a, there is a crossover network size N*(a) at which the 
transition occurs. For TV <C N*(a) the scaling behavior of £ op t is in the strong disorder regime, 
with 4 pt ~ TV 1/3 for ER networks and for SF networks with A > 4, and £ opt ~ jv (A - 3)/(A_1) for 
SF networks with 3 < A < 4. For TV S> TV* (a) the scaling behavior is in the weak disorder regime, 
with £ opt ~ In TV for ER networks and SF networks with A > 3. In order to study the transition we 
propose a measure which indicates how close or far the disordered network is from the limit of strong 
disorder. We propose a scaling ansatz for this measure and demonstrate its validity. We proceed to 
derive the scaling relation between TV* (a) and a. We find that TV* (a) ~ a 3 for ER networks and for 
SF networks with A > 4, and TV* (a) ~ a (A-1 ' /(A_3) for SF networks with 3 < A < 4. 
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I. INTRODUCTION 

The subject of complex networks has been widely ex- 
plored in the past few years in part due to its broad 
range of applications to social, biological and communi- 
cation systems [J H H 0, H . Ina real world network > 
whether it be a communication network or transport net- 
work, the time n taken to traverse a link i may not be the 
same for all the links. In other words, there is a "cost" 
or a "weight" n associated with each link, and the larger 
the weight on a link, the harder it is to traverse this link. 
In such a case, the network is said to be disordered. 

Consider two nodes A and B on such a disordered net- 
work. In general, there will be a large number of paths 
connecting A and B. Among these paths, there is usually 
a single path for which the sum of the costs ^2 t% along 
the path is minimum and this path is called the "opti- 
mal path." The problem of optimal paths on networks is 
of import since the purpose of many real networks is to 
provide an efficient traffic route between its nodes. 

When most of the links on the path contribute to the 
sum, the system is said to be "weakly disordered" (WD). 
In some cases, however, the cost of a single link along the 
path dominates the sum. In this case every path between 
two nodes can be characterized by a value equal to the 
maximum cost along that path, and the path with the 
minimal value of the maximum cost is the optimal path 



between the two nodes. This limit of disorder is called the 
strong disorder (SD) limit ( "ultrametric" limit) and 
we refer to the optimal path in this limit as the min-max 
path. 

The procedure to implement disorder on a network is 
as follows 0, H, El E3- One assigns to each link i of 
the network a random number rj, uniformly distributed 
between and 1. The cost associated with link i is then 



exp(ar,i 



(1) 



where a is the parameter which controls the broadness of 
the distribution of link costs. The parameter a represents 
the strength of disorder. The limit a — * oo is the strong 
disorder limit, since for this case only one link dominates 
the cost of the path. 

There are distinct scaling relationships between the 
length of the average optimal path £ opt and the network 
size (number of nodes) TV depending on whether the net- 
work is strongly or weakly disordered [To| . For strong 
disorder [T^j, i op t ~ N Vovt , where v op t = 1/3 for Erdos- 
Rcnyi (ER) random networks [llj and for scale-free (SF) 
1] networks with A > 4, where A is the exponent charac- 
terizing the power law decay of the degree distribution. 
For SF networks with 3 < A < 4, u opt = (A - 3)/(A - 1). 
For weakly disordered ER networks and for SF networks 



with A > 3, I 



opt 



In TV. Porto et al. 
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optimal path transition from weak to strong disorder for 
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FIG. 1: Schematic representation of the transition in the 
topology of the optimal path with system size TV for a given 
disorder strength a. The solid line shows the optimal path 
at a finite value of a connecting two nodes indicated by the 
filled circles. The portion of the min-max path that is dis- 
tinct from the optimal path is indicated by the dashed line, 
(a) For TV < TV* (a) (i.e. £ oa < t(a)), the optimal path co- 
incides with the min-max path, and we expect the statistics 
of the SD limit, (b) For TV = TV* (a) (i.e. £ ac = f (a)), the 
optimal path starts deviating from the min-max path, (c) For 
TV S> TV* (a) (i.e. loo S> £*{a)), the optimal path has almost 
no links in common with the min-max path, and we expect 
the statistics of the WD limit. 



2-D and 3-D lattices, and found a crossover in the scal- 
ing properties of the optimal path that depends on the 
disorder strength a, as well as the lattice size L. 

Here we show that similar to regular lattices, there ex- 
ists for any finite a, a crossover network size TV* (a) such 
that for TV <C TV* (a), the scaling properties of the op- 
timal path are in the strong disorder regime while for 
TV 3> TV* (a), the network is in the weak disorder regime. 
We evaluate the function TV* (a). The structure of the 
paper is as follows. In Section^ we derive a scaling ap- 
proach for the transition from weak disorder to strong 
disorder of the optimal path. In Section II I II we present 
simulation results which support the scaling Ansatz. Fi- 
nally, in Section ll VI we conclude with an analytic justifi- 
cation for the scaling of the transition. 



II. SCALING APPROACH 

In general, the average optimal path length £ opt (a) in a 
disordered network depends on a as well as on TV. In the 
following we use instead of TV the min-max path length 



ioc which is related to TV as £, 
hence TV can be expressed as a function of £ Q 



£ opt (oo) ~ TW« and 



(2) 



Thus, for finite a, £ op t(a) depends on both a and £oo. 
We expect that there exists a crossover length £*(a), 
corresponding to the crossover network size TV* (a) , such 
that (i) for £00 <C £*(a), the scaling properties of £ op t{a) 
are those of the strong disorder regime, and (ii) for 



^00 3> £* (a) , the scaling properties of £ opt (a) are those of 
the weak disorder regime. In Fig.^ we show a schematic 
representation of the change of the optimal path as the 
network size increases. 

In order to study the transition from strong to weak 
disorder, we introduce a measure which indicates how 
close or far the disordered network is from the limit of 
strong disorder. A natural measure is the ratio 



W{a) 



(3) 



Using the scaling relationships between £ opt (a) and TV in 
both regimes, and £oc ~ N Vo " % (see Section Ij, we get 



£ x ~ N v °* [SD] 
ln4o~kiTV [WD]. 



^opt(a) 

From Eq. and Eq. J3J it follows 
W(a) 



const. [SD] 

toteo/too [WD]. 

We propose the following scaling Ansatz for W(a), 

f 



W(a) = F 



£*{a) J ' 



where 



with 



F(u) 



const. 



u < 1 



\n(u)/u u 3> 1, 



'(a) 



(4) 



(5) 



(6) 



(7) 



(8) 



The dependence of I* {a) on a can be estimated as fol- 
lows. In the strong disorder limit, the costs on the links 
for any path on the network typically differ by at least 
an order of magnitude. This means that for a min-max 
path of £ links (or length £), if we arrange the costs of 
the links in descending order, then two consecutive costs 
typically differ at least by an order of magnitude. If t\ 
and r-i are the random numbers associated with two such 
consecutive links, with r\ > ri, then the ratio of the costs 
on the links is 



— = exp(aAr), 

T2 



(9) 



where Ar = r\ — r%. Thus, in the case of strong disorder 
we must have aAr ^> 1. Consequently the transition to 
weak disorder occurs when all the links become equiv- 
alent in order of magnitude, i.e., when aAr ~ 1. The 
value of Ar depends on the length of the path. If the 
distribution of random numbers on the min-max path is 
uniform, then Ar ~ l/£ for a min-max path of length £. 
The condition for the transition, aAr ~ 1 is satisfied at 
the crossover length £* (a) which implies that 

f(a)~a. (10) 

Therefore, from Eq. ©, W(a) must be a function of 
too/a. 
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III. SIMULATION RESULTS 

Next we describe the details of our numerical simu- 
lations and show that the results agree with our theo- 
retical predictions. To construct an ER network of size 
N with average node degree (k), we start with {k)N/2 
edges and randomly pick a pair of nodes from the total 
possible N(N — l)/2 pairs to connect with each edge. 
The only condition we impose is that there cannot be 
multiple edges between two nodes. 

In order to generate SF networks, we use the Molloy- 
Reed algorithm Each node is assigned a random 

integer k taken from a power-law distribution 



I") 



where ko is the minimal possible number of links that 
a node possesses. Next, we randomly select a node and 
attempt to connect each of its k links with randomly 
selected k nodes that still have free positions for links. 
The disorder in the link costs is then implemented using 
the procedure described in To obtain i^, we use 
the algorithm proposed by Cieplak et al. Q , modified as 
described in Ref. [lJ3- With this modification we reach 
system sizes of N = 2 16 = 65 536. In order to obtain the 
optimal path for a given realization, we use the Dijkstra 
algorithm We calculate the average optimal path 
^opt(o) by taking the average of the optimal paths over 
10 6 pairs of nodes. 

In Fig.|21we show the ratio W(a) for different values of 
a plotted against £ x /£*(a) = £00/ a for ER networks with 
(k) = 4 and for SF networks with A = 3.5. The excel- 
lent data collapse is consistent with the scaling relations 
Eq. ijfjj. Fig. shows the scaled quantities W(a)u = 
£o P t(a)/£*(a) vs. Inu = In^/i* {a)) = ln(4o/«)> for 
both ER networks with (k) — 4 and for SF networks with 
A = 3.5. The curves are linear at large u = £ 00 /£*(a) : 
supporting the validity of the logarithmic term in Eq. {7J) 
for large u. 



IV. DISCUSSION 

We next develop analytic arguments that support 
Eq. (|10|) . These arguments will lead to a clearer picture 
about the nature of the transition of the optimal path 
with disorder strength. 

We begin by making few observations about the min- 
max path. In Fig. 0] we plot the average value of the 
random numbers r n on the min-max path as a function 
of their rank n (1 < n < l^) for ER networks with 
(k) = 4 and for SF networks with A = 3.5. This can be 
done for a min-max path of any length but in order to get 
good statistics we use the most probable min-max path 
length. We call links with r < p c "black" links, and links 
with r > p c "gray" links, following the terminology of 
Ioselevich and Lyubshin 15] where p c is the percolation 
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FIG. 2: Test of Eqs. @ and 0. (a) W(a) plotted as a 
function of £oo /a for different values of a for ER networks with 
(k) — 4. The different symbols represent different a values: 
a = 8(0), a = 16(D), a = 22(o), a = 32(A), a = 45(+), and 
a = 64(*). (b) Same for SF networks with A = 3.5. The 
symbols correspond to the same values of disorder as in (a). 
The insets show W(a) plotted against log(Zoo/a), and indicate 
for loo <S a, W(a) approaches a constant in agreement with 
Eq. Q. 



threshold of the network [T^. We make the following 
observations regarding the min-max path: 

(i) For r n < p c , the values of r n decrease linearly with 
rank n, implying that the values of r for black links 
are uniformly distributed between and p c , consis- 
tent with the results of Ref. [TtJ • This is shown in 
Fig.H 

(ii) The average number of black links, (£b), along the 
min-max path increases linearly with the average 
path length £00. This is shown in Fig.[SjL. 

(iii) The average number of gray links (l g ) along the 
min-max path increases logarithmically with the 
average path length £00 or, equivalently, with the 
network size N . This is shown in Fig. |3Jd. 

The simulation results presented in Fig. [S] pertain to ER 



4 



W(a)u 



W(a)u 




In u 



FIG. 3: (a) Plot of W{a)u = £ op t(a)/t (a) = £ op t(a)/a vs 
In it = ln(£oo/£*(a)) = ]n(£ao/a) for ER networks with (k) — 4 
for different values of a. (b) Plot of W(a)u — £ op t(a)/£*(a,) = 
£opt(a)/a vs Inu — \n(£oo/£*(a,)) = ]n(£oo/a) for SF networks 
with A = 3.5. The values of a represented by the symbols in 
(a) and (b) are the same as in Fig. 



networks; however, we have confirmed that the observa- 
tions (ii) and (iii) also hold for SF networks. Next we 
will discuss our observations using the concept of the op- 
timal spanning tree. The optimal spanning tree (OST) is 
a subset of links of a connected graph which provides an 
optimal path from node A (which serves as the root of 
the tree) to any other node on the graph. When the total 
weight of this path is dominated by the largest weight of 
the links along the path (strong disorder limit), the OST 
does not depend on the root and is determined only by 
the structure of the original graph and a particular real- 
ization of the disorder. In this limit, the OST becomes 
identical to the minimal spanning tree (MST) [H[l|. 
The path on the MST between any two nodes A and 
B, is the optimal path between the nodes in the strong 
disorder limit — i.e, the min-max path. 

To construct the MST, we remove links in the descend- 
ing order of their costs r,-. If removal of a link destroys 
the connectivity of the graph, we restore that link. This 
procedure is continued till there are exactly N — 1 links 



FIG. 4: Dependence on rank n of the average values of the 
random numbers r n along the most probable optimal path for 
(a) ER random networks of two different sizes N = 4096 (□) 
and N = 16384 (o) and, (b) SF random networks. 



remaining. At this point the number of remaining black 
links is 



N(k) Pc 



(12) 



where (k) is the average degree of the original graph and 
p c is given by 0] 



Pc 



(k) 



(k 2 - k) ' 



(13) 



The black links give rise to N c disconnected clusters. 
One of these is a spanning cluster, called the giant com- 
ponent. The N c clusters are linked together into a con- 
nected tree by exactly N c — 1 gray links (see Fig.HjJ. Each 
of the N c clusters is itself a tree, since a random graph 
can be regarded as an infinite dimensional system, and at 
the percolation threshold in an infinite dimensional sys- 
tem the clusters can be regarded as trees. Thus the N c 
clusters containing JVj, black links, together with N c — 1 
gray links form a spanning tree consisting of 2V& + N c — 1 
links. 
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FIG. 5: (a) The average number of links (£b) with random 
number values r < p c on the min-max path plotted as a func- 
tion of its length 1^ for an ER network, showing that (£b) 
grows linearly with l^. (b) The average number of links {l g ) 
with random number values r > p c on the min-max path 
versus IniV for an ER network, showing that (i g ) ~ lniV. 
The inset shows the successive slopes, indicating that in the 
asymptotic limit (£ g ) ~ 1.55 In N. 



Thus the MST provides a min-max path between any 
two points on the graph. Since the MST connects N 
nodes, the number of links on this tree must be equal to 
N - 1, so 



N b + N c = N. 
From Eq. JT2J and Eq. (EJ it follows that 

Nr = N I i _ 



(14) 



(15) 



Therefore N c is proportional to N. 

The path between any two nodes on the MST consists 
of lb black links. Since the black links are the links that 
remain after removing all links with r > p Cl the random 
number values r on the black links are uniformly dis- 
tributed between and p c in agreement with observation 
(i) and Ref. 



FIG. 6: Schematic representation of the structure of the min- 
imal spanning tree, at the percolation threshold, with G be- 
ing the giant component. Inside each cluster, the nodes are 
connected by black links to form a tree. The dotted lines rep- 
resent the gray links which connect the finite clusters to form 
the gray tree. In this example N c — 4 and the number of gray 
links equals N c — 1 = 3. 



Since there are N c clusters which include clusters of 
nodes connected by black links as well as isolated nodes, 
the MST can be described as an effective tree of N c nodes, 
each representing a cluster, and N c — 1 gray links. We 
call this tree the "gr ay t ree" (see Fig. |SJ). This tree is in 



ayti 

fact a scale free tree [l!| with degree exponent X g = 2.5 
for ER networks and for scale for networks with A > 4, 
and \ g = (2A-3)/(A-2) for SF networks with 3 < A < 4 
. If we take two nodes A and B on our original net- 
work, they will most likely lie on two distinct effective 
nodes of the gray tree. The number of gray links encoun- 
tered on the min-max path connecting these two nodes 
will therefore equal the number of links separating the 
effective nodes on the gray tree. Hence the average num- 
ber of gray links (£ g ) encountered on the min-max path 
between an arbitrary pair of nodes on the network is sim- 
ply the average diameter of the gray tree. Our simulation 
results (see Fig. |SJd) indicate that 



InN. 



(16) 



Since (£ g ) ~ ln^oo <C ^oo, the average number of black 
links (£h) on the min-max path scales as £oo in the limit 
of large in agreement with observation (2) as shown 
in Fig-E^. 

Now we will discuss the implications of our findings 
for the crossover from strong to weak disorder. From 
observations (i) and (ii) , it follows that for the portion of 
the path belonging to the giant component, the difference 
between consecutive random values Ar is given by Ar = 
Pc/{(-b)- Hence we can approximate the sum of weights 
of such links by the sum of the geometric progression 



exp(fcaAr) 



k=l 



exp(aAr) [exp(ap c ) 
exp(aAr) — 1 



exp(ar* 



(17) 



6 



where r* 



(l/a)exp(-ap c /<4» ffl. Since (4) 



' p c + - exp 
a 



(18) 



Thus restoring a short-cut link between two nodes on the 
optimal path with p c < r < r* may drastically reduce the 
length of the optimal path. When ap c ^> the proba- 
bility that such a link exists is negligible, but the prob- 
ability becomes significant for ap c ~ £oo- Hence when 
the min-max path is of length w ap c the deviation of 
the optimal path from the min-max path becomes signif- 
icant. The length of the min-max path at which the de- 
viation occurs is precisely the crossover length £* (a) , and 
therefore £*(a) ~ ap c . Note that in the case of SF net- 
works, as A — > 3 + , p c approaches zero and consequently 
I* (a) — » 0. This suggests that for any finite value of dis- 
order strength a, a SF network with A = 3 is in the weak 
disorder regime. 



In summary, for both ER random networks and SF net- 
works we obtain a scaling function for the crossover from 
weak disorder characteristics to strong disorder char- 
acteristics. We show that the crossover occurs when 
the min-max path reaches a crossover length £*(a) and 
£*(a) ~ a. Equivalently, the crossover occurs when the 
network size N reaches a crossover size N*(a), where 



N*(a) 



for ER networks and for SF networks with 



A > 4 and TV* (a) ~ a A ~ 3 for SF networks with 3 < A < 4. 
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